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Abstract. Young energetic pulsars will likely be the largest class of Galactic sources observed by GLAST, with many 
hundreds detected. Many will be unknown as radio pulsars, making pulsation detection dependent on radio and/or x-ray 
observations or on blind periodicity searches of the gamma-rays. Estimates for the number of pulsars GLAST will detect 
in blind searches have ranged from tens to many hundreds. I argue that the number will be near the low end of this range, 
partly due to observations being made in a scanning as opposed to a pointing mode. This paper briefly reviews how blind 
pulsar searches will be conducted using GLAST, what limits these searches, and how the computations and statistics scale 
with various parameters. 
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INTRODUCTION 

Pulsar emission mechanism(s) from radio to gamma-rays are poorly understood after 40 yrs of work. Yet despite this 
fact, pulsars are high-precision tools that probe a variety of topics in both fundamental physics and astrophysics. 
EGRET detected pulsed emission from at least 6, and probably 7—9, young pulsars at energies ^ lOOMeV llj]. In 
addition, several tens of the (primarily Galactic) unidentified EGRET sources are likely pulsars [e.g.|2|,|3l|4j]. With its 
much larger effective area and improved angular resolution, GLAST will almost certainly detect hundreds of pulsars. 
The majority of those pulsars will remain unidentified though (and useless as physics tools), unless pulsations are 
detected through 1 ) "folding" of the events using timing ephemerides for known radio pulsars; 2) searches of associated 
x-ray sources to find radio-quiet Geminga-like pulsars yfl, or very faint radio pulsars (a, 01; or 3) "blind" searches of 
the gamma-ray events [8]. This paper discusses the latter, and perhaps most difficult, option. 

The known gamma-ray pulsars have similar characteristics [for a review, seeQ]]. They are mostly non-variable 
sources with fairly flat energy spectra in the gamma-ray regime, but with spectral cutoffs around or just above 1 GeV. 
The flat energy spectra imply that most of the photons detected by the Large Area Telescope (LAT) will be in the 
~ 100-300 Me V range, where the point spread function has a width £ 1°. The gamma-ray pulse shapes are complex 
with two relatively sharp pulses which provide higher harmonic content in searches and make their detection easier. 

Most of the GLAST mission will be spent in sky survey mode, where the whole sky is scanned every ~3 hourfl 
On average, a point in the sky will be within the LAT field of view ~ 1 /6 of the time. Pointed observations would 
increase this fraction to ~ 1/2 (with Earth occultations preventing higher on-source efficiency). This loss in efficiency 
due to scanning (more specifically, the resulting decrease in the number of source events N s during a particular viewing 
period r v i ew ) will make coherent pulsation searches considerably more difficult and less sensitive. 



BLIND SEARCHES FOR GAMMA-RAY PULSARS 

Blind searches of low count rate event data are typically conducted using either Fourier techniques (i.e. binning events 
into a time series, computing one or more FFTs, and analyzing the resulting amplitude and possibly phase spectra; 
[e-g-SUl) or via brute-force epoch folding (i.e. assembling a pulse profile by determining the pulse phase of each 
event from a trial ephemeris and then computing a probability of non-uniformity for the profile; [e.g. [Tol[Tlll ). Optimal 
sensitivity to pulsations comes from searches that treat all of the data in a "coherent" fashion, meaning that the pulsar's 
rotational phase is accurately tracked over the full observation or analysis duration ?Vi ew , from first event to last. 



For details, see the GLAST mission website: [http : //glast ■ gsf c ■ nasa ■ gov) 




FIGURE 1. Left) The number of independent Fourier "bins" (i.e. l/r v [ ew in Hz) a signal with frequency derivative / (bottom) and 
frequency 2 nd derivative / (top) would drift in an observation of duration r v ; ew . To preserve sharp features in pulse profiles, search 
codes need to account for signal drift to a small fraction (1 — 10%) of a Fourier bin. Right) The 95% confidence-level detectability 
of a blind event-folding search with r v ; ew =20days for a pulsar with a Gaussian pulse profile with FWHM=15% in pulse phase. 
The "sigmas" represent Gaussian significance in standard deviations after accounting for the number of trials searched. The arrow 
(representing a factor of ~3 increase in accumulated events during r v i ew and which can be translated anywhere on the plot) shows 
how the significance improves for pointed observations (which are still affected by Earth occultations) as opposed to scanning 
observations. In this case, a 20-day pointing results in an ~8-<7 detection, whereas 20-days of scanning gives a non-detection. 



The real problem with gamma-ray pulsation searches comes from the fact that the sources have very low count 
rates of l(r 7 -l(r 8 photons(>100MeV) 

cm s , which for the LAT corresponds to roughly 0.2—2 events per 
hour. Therefore, very long observations lasting from weeks to years are required to make significant detections. And 
unfortunately, young pulsars are notoriously badly behaved over such timescales. 

Ideally, we would only need to search over the unknown spin frequency / of a pulsar. However, young pulsars 
rapidly spin-down (i.e. have large frequency derivatives /), they exhibit timing noise which manifests itself in searches 
as significant frequency 2 nd derivatives / (and for the noisiest pulsars, many higher order derivatives as well), and they 
occasionally, on month or year timescales, "glitch" and instantaneously change their observed / and /. Figure 1 shows 
the number of independent frequency bins of width 1 /r v ; ew in Hz, corresponding to an accumulated phase error of 
one rotation, over which a pulsar signal will "drift" during a coherent search for a variety of realistic values of / and 
/. Pulsar spin-down becomes important for energetic young pulsars in a day or two, and so all gamma-ray searches 
must account for an unknown /. Timing noise can become important after a few weeks or months. Since gamma-ray 
pulse profiles contain sharp features, and therefore their pulsed signals contain many harmonics, search codes must 
significantly oversample the ~ 1 /T^ evl or ~ l/^vLw spacings m / an d / to maintain optimal sensitivity. 

A potentially bigger problem, as shown by Chandler et al. jU], are positional errors or uncertainties £, which 
cause apparent changes in the observed spin frequency 8f ~ 10~ 6 £mi- a d/io sin 9 Hz, frequency derivative 8f ~ 2 x 
lCr 13 £ ml - a d/iocos 9 Hz s -1 , and frequency 2 nd derivative 5/~4x 10~ 20 'Emi-ad/io sin0Hzs 2 , where £mrad is the 
position error in milliradians, f\o is the spin frequency in tens of Hz, and 9 is the time varying angle between the 
Earth's velocity vector and the line of sight to the pulsar. For candidate GLAST pulsar sources, £ will be milliradians 
(arcminutes) in size, making the blind detection of even stable Geminga-like pulsars much more difficult. 

The intrinsic f, timing noise, glitch frequency, and position error effects place upper limits on the useful T v { evi for 
coherent pulsation searches to durations as short as ^10 days for the most energetic and active pulsars, or a few 
months for older, slower, and more stable pulsars. Because of the scanning nature of the GLAST mission, there will 
be fewer counts accumulated for a particular source over the specified 7vi ew by a factor of ~3 as compared to a pointed 
observation. This decrease in accumulated counts per unit time greatly impacts coherent search sensitivity and is the 
reason why scanning is not the optimal observing mode for pulsar studies. 



Pulsation Detection and Computational Considerations 



The power P in a periodic gamma-ray signal is P ~ 1 + aNf/Nt, where a is a pulse-shape dependent factor (0.4—0.9 
for the known gamma-ray pulsars)[e.g. 12]. The total number of events N t is the sum of the source (pulsed) events 
N s and background (unpulsed) events Nt,. The probability that noise fluctuations produce a certain search power P is 
exponentially related to the power: Prob(P > P a ) oc exp(— P„) iVtriais, where Aerials is the number of independent trials 
searched. Aerials is proportional to the range of / and / searched, but more importantly, to r v 3 iew . 

In general then, to optimize search sensitivity, we want large N s to maximize P, yet a relatively short 7View 
(appropriate for the limits discussed in the last section) to reduce AViais an d the computational complexity. This rule 
exposes another problem for scanning-mode observations of pulsars: if r v ; ew is increased by a factor of ^3 to recover 
the events "lost" compared to a pointed observation (assuming the pulsations remain well-behaved over the longer 
TviewX a detected signal will have the same power as from the shorter pointed observation, yet the probability that the 
signal is significant will decrease by a factor of ~ 3 3 = 27 due to the larger Aerials searched. 

Computationally, the complexity of epoch-folding searches scales roughly as 7^ (i.e. A^ x AWis), while FFT- 
based techniques fare slightly better and scale as !T v 3 iew log TView- Moderate sized computer clusters can currently handle 
epoch-folding searches of hundreds to thousands of events for viewing periods of up to several weeks (see Figure 1). 
For longer r v i ew or more limited computing resources, one can use evolutionary [14] and/or heirarchical [13] tech- 
niques to greatly reduce the computational burden of coherent searches without greatly impacting the overall sensi- 
tivity, //the average pulsation flux is time invariant. Coherent searches with GLAST with r v ; ew ,$,1—2 months should 
be sensitive to pulsars with photon fluxes (>100MeV) of a few times 10~ 8 to a few times 10~ 7 photons cm~ 2 s _1 for 
scanning-mode observations, or as low as ^2 x 10~ 8 photons cm~ 2 s _1 for pointed observations. 

Finally, incoherent searches are possible where the absolute phase of pulsations is not tracked over the full r v i ew , but 
only over much shorter "windows" r w ; n . Incoherent searches are less sensitive than coherent ones for the same r v i ew , 
but they use many fewer computing resources and enable searches with much longer r v ; ew ll9"l, [l5ll . The basic idea limits 
the duration of the coherently processed intervals r w ; n such that a signal with any reasonable / or / stays within a single 
independent frequency bin (of width 1 /r w ; n Hz) or fraction thereof. These windowed analyses are then combined 
without phase information, but with corrections for / and/or / effects occurring between the intervals, and analyzed. 
Initial simulations show that one might detect relatively stable pulsars with photon fluxes (>100MeV) as low as 
1— 2 x 10~ 8 photons cm~ 2 s _1 with r v i ew ~ 1 yr and r w ; n ~ 1 day. Without phase information, however, the folding of 
all events over such long durations to accumulate high-quality pulse profiles may be difficult or even impossible. 

In summary, instabilities of young pulsars, source positional uncertainties, and the default survey/scanning mode for 
the mission will make blind pulsation searches with GLAST much more difficult than early estimates had predicted. 
It is therefore likely that GLAST will blindly discover only tens of new pulsars rather than hundreds. 

Acknowledgements: Thanks go to Paul Ray, Julie McEnery, Mallory Roberts, Maura McLaughlin, and Kent Wood 
for useful discussions. 
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